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ABSTRACT 

The universality observed in gravitational wave spectra of non-rotating neu- 
tron stars is analyzed here. We show that the universality in the axial oscillation 
mode can be reproduced with a simple stellar model, namely the centrifugal bar- 
rier approximation (CBA), which captures the essence of the Tolman VII model 
of compact stars. Through the establishment of scaled co-ordinate logarithmic 
perturbation theory (SCLPT), we are able to explain and quantitatively predict 
such universal behavior. In addition, quasi-normal modes of individual neutron 
stars characterized by different equations of state can be obtained from those of 
CBA with SCLPT. 

Subject headings: gravitational waves — stars: neutron — stars: oscillations 
(including pulsations) — equation of state — relativity 



1. Introduction 

As is well known, neutron stars of high densities are ideal extra-terrestrial test beds for 
theories of nuclear matters, quark matters and high energy physics (see e.g. Glendenning 
1997, and references therein). For instance, the possible existence of the quark star (see e.g. 
Cheng et al. 1998), a variant of the neutron star, could lend direct support to the theory 
of quark matter (Chodos et al. 1974; Witten 1984; Alcock et al. 1986; Prakash et al. 1990). 
Besides, the relation between the mass and the radius of a neutron star could determine. 
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up to certain accuracy, the equation of state (EOS) of relevant nuclear matter through the 
process of inversion (Lindblom 1992). Therefore, a comprehensive study of observed data 
and statistics of neutron stars is likely to lead to fruitful results in the fields mentioned above. 

Traditional means to gather information about neutron stars have so far relied on elec- 
tromagnetic waves emitted from them. However, with the advent of gravitational wave 
detectors of varioiis designs (see e.g. Hughes 2003; Mason 2004, and references therein), in- 
cluding resonant antennas (e.g. EXPLOPER and NIBOE), ground-based interferometers 
(e.g. LIGO and VIRGO), and the space-based interferometer LISA, it is generally believed 
that neutron stars, as promising gravitational wave emitters, can be observed and analyzed 
in the gravitational wave channel within one or two decades. For example, it has been ar- 
gued that the frequency of detection of gravitational waves emitted in the mergers of binary 
neutron stars could be as high as several hundreds per year in the near future (Belczynski 
et al. 2001; Hughes 2003). Besides, gravitational waves could also be emitted during the 
formation of neutron stars (Lindblom et al. 1998; Fryer et al. 2002). Being spurred on by 
the possibility of inferring the internal structure of neutron stars from the gravitational wave 
signals emitted by them, researchers have been actively examining the peculiarities embed- 
ded in such signals (Andersson & Kokkotas 1996, 1998; Benhar et al. 1999; Kokkotas et al. 



Despite that the evaluation of gravitational waves emitted in violent stellar activities 
such as binary mergers and asymmetric core collapse has indeed posed a grand challenge to 
researchers in numerical relativity, linearized theory of pulsating neutron stars pioneered by 
Thorne & Campolattaro (1967) can still provide useful insight into such complex situations. 
By analogy of radiating electric circuits (Kokkotas & Schutz 1986), hnearized gravitational 
waves are analyzed in terms of quasi-normal modes (QNMs), which are damped harmonic 
pulsations with time dependence exp(ia;i) and are characterized by complex eigenfrequencies 
uj = uj, + \uj; (Press 1971; Leaver 1986a; Ching et al. 1996; Kokkotas & Schmidt 1999; Nollert 
1999). It is believed that the QNM frequencies of a pulsating neutron star can reflect the 
physical characteristics of the star, such as its mass M, radius R and EOS as well (Andersson 
& Kokkotas 1996, 1998; Benhar et al. 1999; Kokkotas et al. 2001). However, Andersson & 
Kokkotas (1998) and Benhar et al. (1999) also noted that the frequency uj-c and the damping 
time r = l/wi of the leading axial and polar w-modes of non-rotating neutron stars follow 
approximately the equations: 
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where Cr, 6r, 0,1, h and ci are model- independent constants determined from curve fitting. 
These universal behaviors seem to undermine the possibility of inferring the characteristics 
of a neutron star from its w-modc QNMs. Instead, they are indicative of a common feature 
shared among neutron stars with different EOS's. 

To seek the physical mechanism underlying the universality in the gravitational wave of 
neutron stars, we showed in a recent paper (Tsui & Leung 2004) that (i) The scaled complex 
eigenfrequencies Co = Muj of w-mode oscillations, to a good approximation, depend only on 
the compactness C = M/R. (ii) Eqs. (1-1) and (1-2) are equivalent to a single formula of a;: 



where a, h and c are complex constants, (iii) The universality indicated by (1-3) is a direct 
reflection of the fact that the the mass distribution inside neutron stars with different EOS's 
can be nicely approximated by the Tolman VII model (TVIIM) proposed by Tolman (1939). 

In this paper wc further pinpoint the crux of the mechanism leading to the universal 
behavior (1-3) for axial w-mode pulsations. Inspired by the success of TVIIM, we propose 
here yet another simple approximation, namely the centrifugal barrier approximation (CBA), 
in which the potential in the axial gravitational wave equation (Chandrasekhar & Ferrari 
1991a) is replaced by a standard centrifugal barrier in the tortoise radius plus a constant 
determined from continuity of the potential at the stellar surface. The CBA potential yields 
the correct asymptotic behavior of the actual potential at the center of the star and is a good 
global approximation as well. The wave function of the axial gravitational wave equation in 
CBA then becomes exactly solvable. It is remarkable that the universal behavior (1-3) can 
still be reproduced with CBA using the tortoise radius of TVIIM star as an input. Hence, 
the observed universality in neutron star axial pulsations is in fact ascribable to (i) the 
centrifugal potential at the center of the star; (ii) the continuity of the potential; and (iii) 
the tortoise radius of TVIIM star, which will be given explicitly in Sect. 2 of our paper. 

To understand quantitatively the universal behavior (1-3), we establish here a scaled 
co-ordinate logarithmic perturbation theory (SCLPT), which is a generalization of the log- 
arithmic perturbation theory (LPT) formulated previously to study QNMs of "dirty black 
holes" (Leung et al. 1997, 1999). Applying SCLPT to CBA and using the tortoise radius of 
TVIIM star, we obtain from first principle the numerical values of a, h and c in (1-3). The 
universal behavior (1-3) is therefore substantiated analytically. 

While the CBA star itself clearly displays universality as mentioned above, we can also 
obtain accurate QNM frequencies of individual realistic neutron stars from those of CBA 
with SCLPT. As a consequence, we are able to consider the EOS-dependence of the scaled 
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QNM frequency Co from the static structure of a star (e.g. mass distribution), thereby 
paving the way for inferring the EOS of a star from its gravitational wave spectra. Thus, the 
universahty and individuahty in QNMs of neutron stars constructed with different EOS's 
are fully explored with SCLPT developed in this paper. 

The organization of our paper is as follows. In Sect. 2, after briefly reviewing the scahng 
behavior of the axial gravitational wave equation, we introduce CBA and show that under 
such approximation QNMs still manifest the generic behavior summarized in (1-3). In Sect. 3 
we formulate SCLPT to study how QNM frequencies are affected by changes in the potential 
and the tortoise radius of the star. In Sect. 4 we apply SCLPT to expand the scaled QNM 
frequency a; as a quadratic function of the compactness of the star and in turn corroborate the 
universal behavior in (1-3). In Sect. 5 SCLPT is applied to investigate the EOS-dependence 
of the scaled QNM frequency Co. We then conclude our paper in Sect. 6 with a discussion 
studying the applicabihty of CBA and SCLPT to wu and trapped modes (see e.g. Kokkotas 
& Schmidt 1999; Chandrasekhar & Ferrari 1991b). Unless otherwise stated, geometrized 
units in which G = c = 1 are adopted in the present paper and numerical results are shown 
for the least damped mode of quadrupole gravitational waves. 



2. Generic behavior of QNMs and CBA 



The eigenvalue equation for QNMs of axial oscillations of neutron stars is given by a 
Regge- Wheeler-type equation (Chandrasekhar & Ferrari 1991a): 

d2 



dr? 



^{u) = 0. 



(2-1) 



Here the tortoise coordinate r* is related to the circumferential radius r by 

r, ^ r e(-''+^y^dr, (2-2) 
Jo 

where e''^^) and e'^^^^ are metric coefficients defined by the hue element ds as follows: 

= -e'^(-)dt2 + e^(")dr2 + r^dO^ + sin^ Odip^). (2-3) 
In addition, the metric coefficient e~'^^^^ is given explicitly by: 

e-^W = 1 - (2-4) 
r 

with m(r) being the mass-energy inside circumferential radius r. The potential V inside the 
star is given by: 



V{r,) = + l)r + Anr^ip - P) - 6m(r)], 



(2-5) 
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where p(r) and P{r) are the mass-energy density and the pressure at a point {r,9,(fi), re- 
spectively. Outside the star, the tortoise radial coordinate reduces to 



= r + 2Mln(^-l)+C. 



(2-6) 



where C is a constant that can be obtained by matching (2-2) with (2-6) at r = i?, and the 
potential V is given by the well known Regge- Wheeler potential (Regge & Wheeler 1957). 



Kw(n) 



1 - 



2M 



6M 



(2-7) 



Eq. (2-1), referred to as neutron star Regge- Wheeler equation (NSRWE) in the following 
discussion, and the outgoing wave boundary condition at spatial infinity together determine 
the QNM frequency of axial w-mode oscillations. We show in Fig. 1 the QNM frequencies 

obtained numerically for neutron stars with different EOS's, including models A (Pandhari- 
pande 1971a) and C (Pandharipande 1971b) proposed by Pandharipande, three models (AU, 
UU and UT) proposed by Wiringa et al. (1988), models APRl and APR2 proposed by Akmal 
et al. (1998), and model GM24 (Glendenning 1997, p. 244). It is clearly shown that both 
the real and imaginary parts of the scaled frequency uj are well approximated by quadratic 
functions in compactness C, as indicated by the dotted lines in the corresponding figures. 
This in turn leads to the universahty discussed in Sect. 1. 

To investigate the physical origin underlying the universality, we have shown that the ax- 
ial mode wave equation displays scaling behavior (Tsui & Leung 2004), and can be rewritten 
as follows: 

" d^ ~ 1 

— +cD2- V(n) 



i^in) = 0, 



where 
and 



l{l + l)f + Airf^ip - P) - 6m{r)^ , f<R/M; 
(1 - 2f-i) [1(1 + l)f-2 - 6f-=^] , f > R/M; 



r — 



m{f) 

P{f) 
p{f) 
v{f) 



r 

M' 

M' 
m{r) 

M ' 
M'^P{r), 

M'p{r), 

iy{r). 



(2-8) 

(2-9) 

(2-10) 
(2-11) 

(2-12) 

(2-13) 
(2-14) 
(2-15) 
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Besides, we also let R = R/M and i?* = R^/M, where i?* = r*(r = R) is the tortoise radius 
of the star. 

Prom (2-8) it is then apparent that the mass-independence of the scaled frequency uj 
follows directly from that of the scaled potential V. Motivated by this and noting that the 
following mass distribution: 

5/r\3 / r \ ^ 



mc{r) = M 



2 \r) ~ 2 \r) 



(2-16) 



conventionally termed as TVIIM in the hterature (Tolman 1939; Lattimer & Prakash 2001), 
is indeed a good approximation of m(r) for neutron stars with various EOSs (see Pig. 2), 
we demonstrated that QNMs of TVIIM (solid circles in Pig. 1) display universal behavior 
observed in realistic stars and interpreted such universality as the consequence of the mass- 
independence of the scaled mass distribution of TVIIM (Tsui & Leung 2004): 



mc(r) 



^c(0 _ 5^3-3 _ 



(2-17) 



In addition to approximating the mass distribution in realistic neutron stars, TVIIM is 
a solvable model whose metric coefficients and can be obtained in closed forms (Tolman 
1939): 



Here ^ = r/R, 



e'^ = (l-5C/3)cos^ 
(f) = (u^i - w)/2 , 



w 



log 



6 + 



3C 



0(^^ — 1) — arctan , 

w{e = 1) , 



c 



3(1 -2C) 



and the pressure is given by: 



1 



C 



\/3Ce-^tan(/)- -(5-3^^) 



AttR^ 

By using (2-18), (2-19) and the definition of r*, we obtain 



R I /(^',C)d^', 
^0 



(2-18) 
(2-19) 

(2-20) 
(2-21) 

(2-22) 

(2-23) 

(2-24) 
(2-25) 
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where 

7(e, C) = { [1 - Ce^(5 - 3^^)] (1 - 5C/3) cos^ 0}"'/' . (2-26) 
Hence, the scaled tortoise radius of the star is 

R.^^ l'HC,C)dC. (2-27) 

Since TVIIM has a mass profile similar to those of realistic neutron stars, we expect dif- 
ferences in R^, between TVIIM and realistic neutron stars would be small. In fact, our 
conjecture is verified by Fig. 3 where R^/M is plotted against M/R for TVIIM and other 
realistic stars. 

As the integral in (2-25) could not be evaluated analytically, we are not able to express 
V{r^) in terms of simple functions of r^, and hence it is not possible to find the analytical 
form of the wave function in NSRWE even inside the star. However, we could expand ^(r*) 
as a power series around the center of the star: 

V{h)^^-^^ + a + 0{fl), (2-28) 

where a is a constant dependent on C and /. It is worthwhile to note that the leading term 
in the expansion of V{r^,) is nothing but a centrifugal barrier in f^,. 

Despite that the leading two terms of the expansion in (2-28) can nicely approximate 
the potential at the center of the star, it usually produces a discontinuity across the stellar 
surface. As QNMs are sensitive to discontinuities, we have to adopt another potential Vc{r*) 
given by 



v>..)^^-^.(.-i)(^-|), (-) 



where i?* is obtained from (2-27). This potential is similar to the one in (2-28) around 
f* = and, in addition, is continuous across the stellar surface. This scheme is termed the 
centrifugal barrier approximation (CBA) in our paper. It is worthwhile to note that under 
CBA the wave function of NSRWE inside the star is simply given by kf^ji{kf^), e.g. for 
/ = 2: 

7, , 3sin(A;r*) 3cos(A;r*) ,r~ ^ ^ , ^ 

i^if*) = ro o - sin(A;r,), < R,. (2-30) 

k''ri kr^ 



Here 

1/2 

(2-31) 



'~2 , ^(^ + 1) A 2^^ /^(^ + i) 6 



k 

Rl V RJ V R^ R^ 
is the scaled effective wave number. 
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In addition to its simplicity, as shown in Fig. 1, CBA also yields QNMs (denoted by 
the solid line in the figure) close to those of TVIIM and realistic neutron stars. Therefore, 
QNMs of CBA manifest the universal behavior (1-3). This remarkable discovery is indeed 
the crux of the present paper. As the scaled potential in (2-29) is completely characterized 
by i? = 1/C and i?*, which is also a function of C, the cause of the observed universality in 
neutron star axial pulsations then becomes obvious and understandable. In a nutshell, the 
universality is attributable to (i) the centrifugal potential at the center of the star; (ii) the 
continuity of the potential; and (iii) the tortoise radius of TVIIM star. 



3. Scaled- Coordinate Logarithmic Perturbation 

After locating the physical origin of the universality in QNMs of neutron stars, we 
further consider the following questions: (i) Can the coefficients a, b and c in (1-3) be 
obtained analytically from CBA? (ii) Can QNMs of reahstic stars be obtained from those of 
CBA? We hold positive views on these challenging issues. In this paper, we aim to develop 
a perturbative study to consider how changes in compactness and EOS could affect the 
frequencies of the QNMs of a star. In fact, the simplicity and generality of CBA directly 
leads to a feasible perturbative analysis for QNMs of neutron stars and in turn provides 
appropriate solutions to these two questions. 

As discussed above, axial oscillations of neutron stars are described by NSRWE with a 
potential term dependent on the distribution of energy-mass density and the pressure inside 
a star. QNMs of these oscillations are the eigen-solution to (2-1) and, in addition, they are 
regular at the origin and satisfy the outgoing boundary condition at spatial infinity. It is 
well known that the wave function of QNMs diverges at spatial infinity and is not amenable 
to standard perturbation theory (sec, e.g. Leung et al. 1999). To this end, we will generalize 
the logarithmic perturbation theory for QNMs, previously formulated by Leung et al. (1997, 
1999) to study how the QNMs of a black hole respond to static perturbations such as a static 
mass shell, to consider QNMs of neutron stars. In this case, the scaled NSRWE (2-8) that 
makes use of the scaled coordinate f is obviously more amenable to perturbative expansion 
since outside the star V(r) depends only f and is independent of the stellar mass. Therefore, 
the method developed here is referred to as the Scaled-Coordinate Logarithmic Perturbation 
Theory (SCLPT). In the subsequent discussion, we will show that shifts in scaled QNM 
frequencies can be expressed in terms of integrals with finite domains of integration. 
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3.1. Pert ur bat ive expansion 

To formulate a pcrturbative expansion for the QNM frequency, we introduce a formal 
perturbation parameter, to measure the departure of the stellar configuration from an 
unperturbed one, which has a scaled circumferential (tortoise) radius Rq {R^o), a scaled 
potential Vo(f*), and a QNM with a scaled frequency cDq. Analogous to other perturbation 
theories, we will assume that the QNM wave function is known. This assumption does 
not pose any problems to our current study because we will use the CBA potential as the 
unperturbed system and, as mentioned above, the wave function can be obtained analytically. 
Besides, R^q and Vo(f*) can be obtained from (2-27) and (2-29), respectively. The subject 
of interest in this paper is to obtain the QNM frequency of a perturbed star whose potential 
function, as a function of f*, can be expressed as: 

= Voir*) + l^Viih) + //V2(r*) + ... , (3-1) 

for f* < -R*(/i). Here -R(/x) and -R*(/i) are the circumferential and tortoise radii of the 
perturbed star, respectively, which are also functions of /i, with R{ii — 0) — Rq and i?*(/x = 
0) = R^Q. 

Outside the star, it is clear from (2-9) that V, as a function of f, is independent of /i 
and hence AV = for f > R{n). It is therefore more convenient to adopt the scaled tortoise 
coordinate f* and the scaled circumferential radius f to describe wave propagation inside 
and outside the star, respectively. Each of these radial coordinates has its own advantage. 
The scaled tortoise coordinate description casts the NSRWE into the standard Klein-Gordon 
equation form, whereas the use of the scaled circumferential radius simplifies the form of V 
outside the star. In our formalism, we will apply these coordinates in different regions of 
interest. 

In the following we seek a power series expansion for the QNM frequency a)(//) of the 
perturbed star: 

a)(//) = tuo + /Ut^i + A*^'^2 + ••• ■ (3-2) 

We will find a general expression for a;„ (n — 1, 2, 3, ...) and explicitly determine the leading 
two expansion coefficients, cui and 0)2. 

3.2. Connection formula 

Inside the star we adopt the scaled tortoise coordinate f* as the independent variable 
and consider the solution to (2-8) as a function of three independent variables f*, a) and /i, 
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namely: 

ip ^'^_{f^,u;,li). (3-3) 

This solution is specified by the regularity condition at the origin and holds for f* < R^. 
Outside the star {f > R), the scaled circumferential radius is used instead and we regard 
the solution to (2-8) there as a function of f and a). In particular, ijj^{r,ijj) is a solution 
that satisfies the outgoing wave boundary condition at spatial infinity. Outside the star, 
NSRWE is identical to the Regge- Wheeler equation of black holes and therefore ilj^{r,uj) 
can be obtained from the Leaver's series solution originally developed for determination of 
QNMs of black holes (Leaver 1986b; Liu 1997; Leung et al. 1999). 

For each /i, we can evaluate the QNM frequency uj by matching the logarithmic deriva- 
tives of ip+ and at r = R{ii) = Rq + iiRi + /x^i?2 + • • •, i-e. 

/_(^.(//),a;(//),/i) = /+(%), a;(/i)), (3-4) 

where 

f-{f^,Cj,ii) = ^-J-: — (3-5) 



1 dilj+ 

(3-6) 



(r-2) djj^ ^ 

It is interesting to note that f± are matched at the scaled radius i?(/x) (or i?*(/x)) that would 
vary with /j,. 



3.3. Expansion of logarithmic derivative 

To generate an expansion for the QNM frequency u}{iJ,), we firstly expand /_ as a power 
series in the perturbation parameter /x: 

f_{f*,Lb{ix),ix) = h{r*)+ nh{f*)+ 1? h{f + , (3-7) 

where f_{r^,UQ,^ = 0) = fo{r*) and we have suppressed the dependence on ujq in the 
expansion coefficients fi{r*) {i = 1,2,3,...). By transforming (2-8) into the form of the 
Riccati equation: 

f_ + f_+u^-V{h) = 0, (3-8) 
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making use of the expansions (3-1) and (3-2), and comparing equal powers of /i, we show 
that 

f:, + 2fofn + 2UoC^n^Un: ^ = 1, 2, ... . (3-9) 

Here C/i(f*) = Vi(r*) and for n > 2, 

n-l 

Un{r*) = K(^*) - ^[Mr*)fn-i{r*) + ojiOJn-i]: (3-10) 

i=l 

and we have used the prime to symbohze differentiation with respect to f*. Eq. (3-9) can 
be solved by introducing the integrating factor exp[2 J^* dyfo{y)] — V'o(^*)' resulting in an 
integral expression for fniR*o)' 

fn{R*o)tpl{R*o) ^ dh[Un{h) -2uoUn\i>o{r*)- (3-11) 

Jo 



Secondly, as the logarithmic derivatives are to be matched at f* = i?*(/x) = R^o + f^R*i + 
A*^-R*2 • •, we expand f-{R^{iJ,), 00(11), n) and /+(i?(//), ci;(//)) in power series of /i, yielding: 



f_{R^{li),ij{li),li) 



and 



MMi^)) + f^MRM) + i^'f2{R*{fi)) + 

dfn /i^ d^/o 



/o(-R*o) + /i 



+ 



+ 



/i(-R*o) 



d/i ' 2! d/i2 



/2(-R*o) +1^ 



dfi ' 2! d/x2 
d/2 , /x^dV: 



+ 



d// 2! d//2 



+ ... 



+ ... 



+ 



(3-12) 



U{R{pi),u;{pi)) = U{R{f,),u;{pi)) 
= f+{Ro,u}o) + 1^ 



di+ w^d^f 
d/i 2! d//2 



+ 



+ ... . 



(3-13) 



Here it is understood that all derivatives with respective to /i are evaluated at /i — 0. Hence, 
by comparing equal powers of /i in (3-12) and (3-13), fn{R*o) can be expressed in terms of 
derivatives of fi {i < n) and /+ as follows: 



fn{R*o) — 



1 d"/+ 



1 d"/o 
n! d/i'^ 

On, 



+ 



/l 



— 1)! d/i 



n-l 



H h 



d/x 



(3-14) 
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where, by definition, ©„ is equal to the sum in the square bracket. Furthermore, by intro- 
ducing another expression A„ defined by: 

U{Ro + iiRi + --- + n''-^Rn^i,Cuo + i2Cui + --- + ii^'-^u^^i) - f+{RQ,CjA (3-15) 

J n 

where [F(/x)]„ indicates the nth-order term of a function F(/x), it is readily shown that 

Id"/" ( d d\ 

+^ /^(cDo,i?o). (3-16) 



(3-17) 



As a result, we obtain a formal expansion for the scaled frequency a), which reads: 

OJji — — , 

2a;o(^o|^o) 
with 

(HUnli^o) = 2 *° dnUn{f*)'4^oir*) + (^©n - A„ - Rn^^ ^l{R*o). (3-18) 



and 



(^ol^o) = r i^lirMr. + ^^i^- (3-19) 



Despite the simplicity of (3-17), which is akin to standard perturbation formulas in 
quantum mechanics with playing the role of the norm squared of a quantum state, 

the emergence of terms like f/„, A„ and 0„ reveals the achievement underlying this formula. 
It is worthy of remark that in (3-18) there are three different contributions to the frequency 
shift ujm namely an integral over the interior of the star and two surface terms originating 
from /+ and /_, respectively. In the following discussion, we will work out explicit expressions 
for the first and second order results. 



3.4. First and second order frequency shifts 

For the case n = 1, the expression of t/i(f*) is trivial, while 0„ and A„ are given by: 

e, = ^ = (3-20) 
dfi dr* d/x 

Ai = 0. (3-21) 

As a result, it is clear that 

iHUSo) = ^ dnVi{n)i>Uf.) + {r*i^ - ^i^) V^o(^*o), (3-22) 
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and hence cDi can be readily obtained from (3-17) and (3-19). 
For the case n = 2, it is straightforward to show that 



IdVo d/i 

2 dfj,'^ d/i 

d/o^ , ldVo^2 , dA~ . 

dr=^ 2 dr^ dr^, 

- 2^ia^ + ^^'^^afa^ + 2^^^- 

Direct substitution of these results into (3-18) leads to the second-order term: 

(^0|^2|^0) 

V2{f^) - /i (f*) - ojf 'tpQ{f^)df^ 



df* * 2 df^ df* 



+V^o(^*o) 
-^o(^*o) 



--R*o 



2 acD2 



orouj 



2 af2 



(?f 



f=Ro 



and the second-order frequency change follows directly from (3-17). 



(3-23) 



(3-24) 



(3-25) 



(3-26) 



After succeeding in deriving the first and second order shifts in the eigenfrequency of 
QNMs of neutron stars, we will apply relevant formulas to answer the questions posed in the 
beginning of this section. 



4. Universality in QNMs 

Being a good global approximation to realistic stars, CBA also demonstrates the uni- 
versal behavior summarized by (1-3). In fact, the solid line in Fig. 1, representing the QNM 
frequencies of CBA, is close to the best quadratic fit to those of realistic stars (the dotted 
line). Therefore, we expect that the universal behavior displayed by realistic neutron stars 
can be understood from the QNMs of CBA and Eq. (1-3) can be deduced from CBA as well. 
Motivated by this conjecture, we evaluate the QNM frequency of CBA with SCLPT. 
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In SCLPT, the QNM frequencies for CBA stars with different compactness can be 
obtained by considering compactness C as the perturbation parameter. In this case, the 
formal expansion parameter /x = C — Cq, where Cq is the compactness of a reference CBA 
star whose QNMs are known. Hence, to second order in C — Cq, u;{C) is approximately given 
by: 

u{fi) =Coo + iC- Co)ui + {C- CofCj2 , (4-1) 

where ujq is the QNM frequency of the reference CBA star. The first and second order 
shifts in QNM frequencies are proportional to C — Co and (C — Cq)^, respectively. Under this 
approximation, the QNM frequency is expressible in terms of a quadratic function of C: 



with 



a) = aC^ + 6C + c , (4-2) 

a = a'2; (4-3) 
h = ui- 2C0O2; (4-4) 

C = LVq- CqLOi + Cli02. (4-5) 

It is obvious that Eq. (4-2) is in prefect agreement with the universal behavior (1-3) dis- 
covered numerically by Andersson & Kokkotas (1998); Benhar et al. (1999); Tsui & Leung 
(2004). 

To gauge the accuracy of the second-order SCLPT mentioned above, we apply it to a 
reference CBA star with compactness M/R = 0.2. As shown in Fig. 4, the results obtained 
from SCLPT (represented by the solid line) arc good approximation of the exact QNMs 
(represented by the stars). This clearly demonstrates the validity of SCLPT. In addition, the 
pcrturbative results also faithfully demonstrate the universal behavior displayed by realistic 
stars (represented by the dotted line), and the numerical values of a, b and c obtained from 
(4-3), (4-3) and (4-3) respectively are in nice agreement with the those obtained from the 
best quadratic fit to the QNMs of the realistic stars (see Table 1 for reference) . Hence, the 
universality in axial pulsations of neutron stars is fully understood and predicted analytically. 

The technical details of the perturbation scheme yielding uji, u;2 and hence the constants 
a, b and c are as follows. To obtain the first and second order shifts, we have to evaluate the 
all quantities appearing in (3-22) and (3-26). Specifically, ^1 and R2 are given by: 

(4-6) 

(4-7) 
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Analogously, R^i and i?*2 can be obtained: 



dC 



C=Co 



1 r 



R*2 — — 



1 

^0 Jo 
1 / d^^. 



dCo Co 



(4-8) 



2 I dC2 



C=Co 



2Co Jo 



Co aco 



+ 



*^0 



(4-9) 



On the other hand, as Vc contains explicit dependence on i?* and R, V\ and Vi can be found: 



Vl(r*) 



dC 

' C=Co 

2/(/ + l)7?,i 

^^0 



2l{l + l)Co - 3(/^ + / + ^)Cl + l2Cl 



C=C(, 



6/(/ + l)i?2^ 4/(/ + 1)R,2 



+ 



i?2o 



+ 2l{l + 1) - 6(^2 + i + 3)Co + 36Co', (4-10) 



which are constants independent of f*. 



It is easy to see that T4(^*) {t^ — 1,2,3,.. .) are indeed all f*-independent under CBA. 
This fact greatly simplifies the perturbation calculation in CBA and renders the integrals 
involved in the evaluation of cDi and 0)2 exactly solvable. For example, noting that "00 (^*) — 
kr^ji{kf^), we have 



R*o 



and 



dr*yi(r*)V'o(r*) 



ViR. 



*° W [jfi^) + Jf-M - (2^ + ^)^Jii^)ji-^i^)}M. (4-11) 



^o(^*) JO 
{Vi - 2u}ooJi)R*o 



Jliix) (2/ + l)j;_i(a;) 



1 + 



xji{x) 



(4-12) 



x=kR*o 
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Besides, we have also apply SCLPT to evaluate QNMs of other realistic stars. As an 
example, we show here the result for a star constructed with APR2 EOS. As demonstrated 
in Fig. 5, the second order result of SCLPT is again a good approximation to the exact 
numerical results. This provides independent corroboration to the validity of SCLPT. 



5. Individuality of QNMs 

After establishing the universality in axial pulsations of neutron stars, we turn our 
attention to the individuality of such QNMs. As clearly shown in Figs. 1 and 6, QNM 
frequencies of individual realistic neutron stars in general deviate slightly from the universal 
curve (1-3). In the present paper wc aim to evaluate QNMs of realistic stars from those 
of CBA star with SCLPT developed here. For a realistic star whose compactness C and 
potential V"(f*) are known, we compare it with a CBA star with the same compactness C. 
Therefore, we have Vq = ^ = K + ^i; and i?* = R^q + where R^q is the tortoise 
radius of the CBA star. However, in this case the scaled circumferential radii of the two 

stars in consideration are the same. In other words, i?„ = for n = 1,2, Besides, it is 

convenient to take the formal expansion parameter /i — 1 in the current situation. As both 
= and R„i = for n = 2, 3 . . ., the matrix elements ('0o|t^i|'0o) a-nd ('0o|t^2|'0o) can be 
simplified and respectively take the form: 

iHUilA) = r" dhVi{h)ip'oir*) + R^i^ipliR^o), (5-1) 
Jo dr* 



(^0|^2|^0) 

" [fi{r*)+^!] 4'(r*)df. 



+^^0(^*0) 



2 df2 * df* 



f*=R*Q 



f=Ro 



(5-2) 



The first and second order shifts in QNM frequency then follow directly from (3-17), (3-18) 
and (3-19). 

Figure 6 shows QNMs obtained from the perturbation scheme outlined above for stars 
constructed from various EOSs with a common compactness C — 0.20. The unfilled, dark, 
and grey symbols respectively indicate the exact, the first and the second order results, 
while the star is the unperturbed frequency of CBA. It is clearly shown that the second 
order perturbation results can nicely approximate the exact ones. We have also applied 
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SCLPT to stable neutron stars with a larger (or smaller) compactness and found that the 
second-order results are indeed reliable. In fact, the first order formula readily suffices to 
yield accurate prediction of the shift in Reci;. 

On the other hand, we have to point out the limitations of SCLPT. As in other cases, 
perturbation theory is likely to break down in the presence of instabilities. Therefore, we 
expect that the accuracy of perturbativc results obtained from SCLPT gets worse when the 
star becomes unstable against perturbation. We verify this point by applying SCLPT to 
neutron stars constructed with GM24 EOS (Glendenning 1997, p. 244), which are known 
to be unstable if the compactness is greater than 0.21. As clearly demonstrated in Fig. 7, 
where QNMs of CBA and GM24 stars with C = 0.17,0.19,0.20 are shown, the deviation of 
the second-order result from the exact value increases as the compactness approaches the 
border of stability. However, as far as Re a; is concerned, both the first and second order 
results are satisfactory. 



6. Conclusion and discussion 

The main achievement in the present paper is the discovery that the universality in the 
QNM frequency of axial pulsations of neutron stars in fact originates from the CBA, under 
which the potential term in the scaled NSRWE is essentially a centrifugal barrier. There 
are two parameters in the scaled potential, namely the scaled circumferential radius R and 
tortoise radius i?* of the star. While the former is just the inverse of the compactness C, the 
latter can be determined from that of TVIIM. These two parameters completely determine 
the QNMs of CBA and in turn give rise to the observed universality. 

In order to consider how the the physical characteristics of a neutron star could affect the 
frequencies of its QNMs, we have also developed a systematic perturbativc scheme SCLPT, 
which is designed to cope with the divergence in QNM wavcfunction at spatial infinity. 
Direct application of SCLPT to CBA then successfully predicts the universality in the QNM 
frequency of realistic neutron stars. The advantage of application of the scaled coordinates 
(f and f*) in our study becomes manifest in this regard because the potential V outside 
a star is independent of its detailed internal structure save its compactness, leading to the 
observed universahty (Andersson & Kokkotas 1998; Benhar et al. 1999; Tsui & Leung 2004). 

On the other hand, SCLPT can also predict small deviations in QNM frequencies from 
those of CBA for individual realistic neutron stars. This opens possibilities for researchers to 
fink the QNM frequencies observed from a distant star with its internal structure and EOS 
as well. We are currently working on the inverse problem of SCLPT, namely inferring the 
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internal structure of a pulsating neutron star from its gravitational wave spectrum. Relevant 
progress in this direction will be reported elsewhere in due course. 

So far we have used the least-damped w-mode to illustrate the principle and accuracy 
of our method. However, the results mentioned in the present paper are general and hold for 
modes of higher orders (i.e. modes with frequencies higher than that of the least-damped 
one). More interestingly, we note that the wn-mode (sec e.g. Kokkotas & Schmidt 1999, and 
references therein) also displays similar universality. As shown in Fig. 8, where the real and 
imaginary parts of the scaled QNM frequencies of the uin-niodes of various realistic stellar 
models, TVIIM and CBA are plotted against the compactness, the universality summarized 
in (1-3) is clearly exhibited. It is remarkable that the corresponding values obtained from 
application of second-order SCLPT to CBA (the solid line) successfully capture the essence 
of such universal behavior. To further examine the validity of CBA and SCLPT in this case, 
we show in Fig. 9 the scaled QNM frequency of a u^n-mode of a CBA star (the star symbol) 
and other realistic stars (the unfilled symbols) with a common compactness C = 0.2, the 
first and second order results obtained from applying SCLPT to the CBA star (respectively 
denoted by the corresponding dark and grey symbols). It is obvious that SCLPT can indeed 
reproduce accurate eigenfrequencies for the wn-mode. 

Despite that we have demonstrated here the validity of CBA and SCLPT for ordinary 
tf-modes and wji modes; we have to caution that the trapped mode is an exception to 
our method. As is well known, trapped modes usually exist in ultra-compact stars with 
compactness greater than 1/3 and have small decaying rates (Chandrasekhar & Ferrari 
1991b). The physical origin of such modes is the development of a local minimum in the 
potential V{r^) inside the star when C > 1/3. Hence, gravitational waves are trapped in the 
potential minimum and acquire much longer lifetime. The analysis proposed in the present 
paper relies on the vahdity of the CBA model whose potential term V{r^) is obviously a 
monotonic function inside the star. In fact, the CBA potential deviates significantly from 
those of TVIIM (or other ultra-compact stars). As shown in Fig. 10, trapped modes of 
CBA are no longer good approximation to those of TVIIM if C > 1/3 and therefore it is 
not possible to obtain the trapped modes of TVIIM (or other ultra-compact stars) from 
application of SCLPT to CBA. We conclude our paper with this remark. 

We thank J Wu for discussions. We also express our gratitude to an anonymous referee 
for drawing our attention to the Wu and the trapped modes. Our work is supported in part 
by the Hong Kong Research Grants Council (grant No: CUHK4282/00P and 401905) and a 
direct grant (Project ID: 2060260) from the Chinese University of Hong Kong. 



- 19- 
REFERENCES 

Akmal, A., Pandharipande, V. R., & Ravenhall, D. G. 1998, Phys. Rev. C, 58, 1804 
Alcock, C, Farhi, C. E., & Olinto, A. 1986, ApJ, 310, 261 
Andersson, N., & Kokkotas, K. D. 1996, Phys. Rev. Lett., 77, 20 
— . 1998, MNRAS, 299, 1059 

Belczynski, K., Kalogera, V., & Bulik, T. 2001, ApJ, 572, 407 
Benhar, O., Berti, E., & Ferrari, V. 1999, MNRAS, 310, 797 
Chandrasekhar, S., & Ferrari, V. 1991a, Proc. R. Soc. A, 432, 247 
— . 1991b, Proc. R. Soc. A, 434, 449 

Cheng, K. S., Dai, Z. G., Wei, D. M., & Lu, T. 1998, Science, 280, 407 

Ching, E. S. C., Leung, P. T., Sucn, W. M., & Young, K. 1996, Phys. Rev. D, 54, 3778 

Chodos, A., Jaffe, R. L., Johnson, K., Thorne, C. B., & Weisskopf, V. F. 1974, Phys. Rev. 
D, 9, 3471 

Fryer, C. L., Holz, D. E., & Hughes, S. A. 2002, ApJ, 565, 430 

Glendenning, N. K. 1997, Compact Stars - Nuclear Physics, Particle Physics, and General 
Relativity (Springer, NY) 

Hughes, S. 2003, Ann. Phys., 303, 142 

Kokkotas, K. D., Apostolatos, T. A., & Andersson, N. 2001, MNRAS, 320, 307 
Kokkotas, K. D., & Schmidt, B. G. 1999, Living Rev. Rel., 2, 2 
Kokkotas, K. D., & Schutz, B. F. 1986, Gen. Relativ. Gravitation, 18, 913 
Lattimer, J. M., & Prakash, M. 2001, ApJ, 550, 426 
Leaver, E. W. 1986a, Phys. Rev. D, 34, 384 
— . 1986b, J. Math. Phys., 27, 1238 

Leung, P. T., Liu, Y. T., Suen, W. M., Tam, C. Y., & Young, K. 1997, Phys. Rev. Lett., 78, 
2894 



-20- 



— . 1999, Phys. Rev. D, 59, 044034 
Lindblom, L. 1992, ApJ, 398, 56 

Lindblom, L., Owen, B. J., & Morsink, S. M. 1998, Phys. Rev. Lett., 80, 4843 

Liu, Y. T. 1997, MPhil thesis. The Chinese University of Hong Kong 

Mason, J. W., ed. 2004, Astrophysics Update, Vol. I (Springer-Praxis), 281-310 

NoUert, H.-P. 1999, Class. Quantum Grav., 16, R159 

Pandharipande, V. 1971a, Nucl. Phys A, 174, 641 

— . 1971b, Nucl. Phys A, 178, 123 

Prakash, M., Baron, E., & Prakash, M. 1990, Phys. Lett. B, 243, 175 
Press, W. H. 1971, ApJ, 170, L105 

Rcggc, T., & Whcclcr, J. A. 1957, Phys. Rev., 108, 1063 

Thorne, K. S., & Campolattaro, A. 1967, ApJ, 149, 591 

Tolman, R. C. 1939, Phys. Rev., 55, 364 

Tsui, L. K., & Leung, P. T. 2004, MNRAS, 357, 1029 

Wiringa, R. B., Fiks, V., & Fabrocini, A. 1988, Phys. Rev. C, 38, 1010 

Witten, E. 1984, Phys. Rev. D, 30, 272 



This preprint was prepared with the A AS IM^jX macros v5.2. 



-21 - 




Fig. 1. — The real and imaginary parts of the scaled QNM frequencies of various realistic 
stellar models (including APRl, APR2, AU, GM24, Models A and B, UT and UU), TVIIM 
and CBA (solid line) are plotted against the compactness. The dotted line is the best 
quadratic fit to those of the realistic stars 
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Fig. 2. — The mass distributions of various realistic neutron stars and TVIIM with a common 
compactness C — 0.2 are plotted against r/R. 




Fig. 3. — The figure shows the relationship between R*/M and the compactness for TVIIM 
(solid line) and other realistic neutron stars. 
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Fig. 4. — The scaled QNM frequency of the CBA star is plotted against the compactness 
C. The stars and the sohd hne represent the exact values and the ones obtained from 
applying the second-order SCLPT to a CBA star with C — 0.2, respectively. For purpose of 
comparison, we also show the universal curve (the dotted line), i.e. the best quadratic fit to 
QNM frequencies of the realistic stars considered in this paper. 
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Fig. 5. — The scaled QNM frequency of the APR2 star is plotted against the compactness 
C. The stars and the sohd line are the exact values and the ones obtained by applying the 
second-order SCLPT to an APR2 star with C — 0.2, respectively. 
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Fig. 6. — The star symbol shows the scaled QNM frequency of a CBA star with compactness 
C = 0.2. By applying SCLPT to the CBA star, the first and second order results, respectively 
denoted by dark and grey symbols, for QNM frequencies of realistic stars with the same 
compactness are obtained. For comparison, exact numerical QNM frequencies of realistic 
stars are indicated by the corresponding unfilled symbols. 
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Fig. 7. — The figure shows the scaled QNM frequencies of GM24 stars with C = 
0.17,0.19,0.20. The star symbol indicates the QNM frequencies of the corresponding CBA 
stars. The first and second order results obtained from SCLPT are denoted by dark and 
grey symbols, respectively. For comparison exact numerical QNM frequencies of GM24 stars 
are indicated by the corresponding unfilled symbols. 
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Fig. 8. — The real and imaginary parts of the scaled QNM frequencies of a -uJn-mode of 
various realistic stellar models, TVIIM and CBA are plotted against the compactness. The 
solid line shows the corresponding values obtained from the second-order SCLPT. 
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Fig. 9. — The star symbol shows the scaled QNM frequency of a wn-mode of a CBA star 
with compactness C = 0.2. By applying SCLPT to the CBA star, the first and second order 
results, respectively denoted by dark and grey symbols, for QNM frequencies of realistic stars 
with the same compactness are obtained. For comparison, exact numerical QNM frequencies 
of realistic stars are indicated by the corresponding unfilled symbols. 
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Fig. 10. — The real and imaginary parts of the scaled QNM frequencies of a trapped- mode 
of TVIIM and CBA are plotted against the compactness. 
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Data 


a 


h 


c 


Best fit to realistic stars 
Perturbative result of CBA 


-3.9-i5.6 
-4.63 - 14.86 


2.8 + 11.6 
3.15 + 11.30 


-0.03 + 10.125 
-0.063 + 10.155 



Table 1: The complex coefficients a, h and c in (1-3) are shown. In the first row the values 
are obtained from the best quadratic fit to the QNMs of the realistic stars considered in 
Fig. 1, while in the second row those obtained from the second order result of SCLPT for a 
reference CBA star with compactness M/R — 0.2 are presented. 



